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ABSTRACT 


We discuss a fundamental radar imaging problem for narrowband electromagnetic waves that 
extends the recent results originally obtained in [1, 2] in the scalar, or acoustic, setting. In particular we 
demonstrate the ability to efficiently image three-dimensional convex conducting bodies by using the 
knowledge of the scattered electric field for one fixed monochromatic illumination of the target. In this 
problem our knowledge of the scattered electric wave is understood to be in the form of measurements of 
the amplitude and phase of the tangential components of the radiated electric field (generated by the radar 
target) on a finite two-dimensional array. We also provide a more general version of this streamlined result 
that describes imaging of multiply connected nonconvex objects with the same measurements. 
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1. INTRODUCTION 


We begin with a few introductory remarks that outline what we mean by narrowband 
electromagnetic imaging. The problem concerns estimating the shape and size of an object on the basis of 
known information in the form of measurements of reflected waves taken at some distance from the 
unknown object. Specifically, how do we form a three-dimensional image of a potentially very distant 
object using something like a radar or collection of such sensors? Figure 1 depicts such a measurement 
scenario. 



Figure 1. Data collection scenario. 


The problem is simple: image the scatterer, or target, with a single incident wave with the aid of 
measurements of the scattered wave taken on a finite two-dimensional array. For example, the scattered 
wave is a frequency-dependent, i.e., narrowband, induced voltage response created by the illumination of 
the target. Although the underlying physics and mathematics that lead us to the imaging algorithm can 
appear complicated, the algorithm we will ultimately employ is extremely easy to understand and 
implement. 

The imaging algorithm may be summarized as follows. The measured electric field, as a function of 
two variables on the array and for a chosen origin of the problem, is decomposed into its generalized 
Fourier modes, and a combined one-dimensional power spectrum of this signal is generated. The critical 
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point in the power spectrum where the combined spectral energy transitions from a highly oscillatory 
behavior and begins a super exponential decay is then found. This index is determined by computing the / 2 
norm 1 of the accumulated energy of these modes and defining the cutoff at ninety percent of the total 
spectral energy. The index cutoff is converted to a length through division by the wave number (2it/X) of 
the narrowband waveform used in probing the target. This length represents the radius of the smallest 
sphere, centered at the origin we previously specified, which must contain the unknown scatterer. The 
process is then repeated with new center, or origin, and terminated when enough centers of interest have 
been considered. Finally, the intersection of all the computed spheres is computed and the resulting 
common region is taken as an approximation to the convex scatterer. 

We mention that the power spectrum possesses three very pronounced and particular behavior 
regimes as a function of increasing mode number. Figure 2 illustrates the three regions of interest and 
shows the combined power spectra, labeled as | (3 n | for a candidate example. The first regimes takes place 
within 0 < n < 45, the second roughly within 45 < n < 55, and the final region for all n beyond this. 


1 . 


i.e., L 2 


norm of a discrete collection of complex numbers. 
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Figure 2. Oscillation and decay of the combined power spectra. 


First, there is a region of highly oscillatory modes, followed by a region of extremely rapid decay, 
finally followed by a region of essentially null values. The index at which this power spectrum begins this 
rapid decline defines the radius of the smallest sphere with a specified origin, which must contain the 
scattering object. 2 By considering several origins we construct several such spheres, each of which must 
contain the scatterer. Hence, by considering the intersection of all such spherical regions we arrive at an 
approximation to a three-dimensional image of the assumed target. Figure 3 provides an illustration of this 
process. 


2. For the moment, we mention that the reason for this is intimately related to the unique analytic continuation of the 
measured electric field on an array. We will revisit this in the section treating the supporting theory behind our main 
result and defer further discussion until such time, due to its technical nature. 
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Figure 3. An example of two-dimensional projection of the intersection of three spherical-containing regions and the 
true convex target. 


The remainder of this report reviews the physics that govern this problem and how that physics and 
some mathematical analysis lead to a well-understood and fully justified target characterization algorithm. 
Moreover, this treatment suggests additional innovative object characterization techniques through its 
description of the fundamental scattering processes at work. These results will be presented in forthcoming 
articles. 
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2. THE PHYSICAL MODEL 


We begin with a brief overview of the fundamental inverse scattering problem at the heart of our 
application. Our imaging technique will employ novel processing of measurements of electromagnetic 
waves. Hence, to illustrate the origins, success, and limitations of the algorithm, we begin with a brief 
review and characterization of what we term the forward, or direct, scattering problem: given a well- 
defined scattering object and a specified incident wave, we formulate the response to that incident 
illumination and analytically predict how such a wave would interact with the scatterer. 

We start with the time-dependent Maxwell’s equations in free space, i.e., in a subdomain of £2 c R 
in which there are no electric or magnetic current densities, for which the electric field £ and the magnetic 
field B satisfy 


V a £(x, t) + p 0 3 ( B(;t, r) = 0 , 


( 1 ) 


V a B(x, t) — e 0 d t E(x, t) = 0. 


( 2 ) 


Here, x = (x x , x 2 , * 3 ) e R } is our spatial variable and (e R denotes time. The constant quantities 
p 0 and e 0 are, respectively, the magnetic permeability and electric permittivity of free space. We briefly 
mention that the electric and magnetic fields are complex-valued vector quantities as functions of both 
position and time, so that £, B: R 3 x R— »C 3 . Additionally, V a is the standard curl operation and d , 
denotes differentiation in time. 

Rather than dealing directly with the problem in the time domain, we will recast this system of six 
coupled partial differential equations into one that depends on frequency. On the basis of this system, we 
arrive at our narrowband, or fixed-frequency, imaging result. 

To accomplish this goal, we employ the Fourier transform mapping functions, or distributions, in 
time into ones depending on angular frequency w. In this way we define frequency-dependent electric and 
magnetic fields as 


E(x,w) = je lwl £(x,t)dt and H{x,w) - Je B(x, t)dt , 
R R 

whose inversions are given by 


( 3 ) 
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( 4 ) 


E(x, t) = — \e E(x,w)dw and B{x,t) = — [e H(x,w)dw. 

IK J 2n J 

R R 

We note that i = represents the imaginary unit. 

Employing the Fourier transform essentially converts differentiation in time into multiplication in 
frequency. Hence, we arrive at the frequency-dependent Maxwell system 


V a E(x, w) - iw\i 0 H(x, w) = 0, (5) 

V a H(x, w ) + iw£. Q E(x, w) = 0 . (6) 

We now take a moment to distinguish between two fundamentally different types of electric and 
magnetic fields. By this we mean outgoing, or scattered or radiating waves {Ef, H s ), and incoming, or local 

incident plane waves ( E ‘, H l ). We define the total fields we would measure with any physical device as the 
sum of these two types of waves. Namely, we will use, and so we define here, the total fields E and H as 


E(x) = E\x) + E\x), H(x) = H\x) + H\x). (7) 

Any outgoing pair of electric and magnetic fields must satisfy the so-called Silver-Miiller radiation 
condition, which we recall is given by either 


lim (H s ax - rE s ) = 0, r = |j 


( 8 ) 


or 


lim ( E s a x + rhf) = 0, r = |jc| 


(9) 
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Equivalently, E? = (E\, E? 2 , Ef 3 ) and H s = (Hi, H 2 , H 3 ) means that each Cartesian coordinate 

v 3 s 3 

component Ey.R —> C and Hy.R = 1,2,3, of the two vector fields E? and H s satisfies the 

Sommerfeld radiation condition, which is 


lim ( d r -ik)E = lim ( d r ik)H = 0, r = |jc| . (10) 

r—>00 J r —> co J 


The Silver-Miiller and Sommerfeld radiation conditions imply the asymptotic behaviors 


ik\x\ 

e \ x ) = ^(^w+ocur 1 )), w->oo, (id 

and 

«*W 

= -^-(H„Cx) + 0(\4 "')), W-*~, (12) 

where x = (x/|x|) and |x|“ = x] + x 2 + x\ . Here, the quantities E^ and denote the far-field patterns, 
or scattering amplitudes, of the outgoing waves E s (x) and H s (x). We emphasize that for most remote 
electromagnetic sensing experiments, these are the quantities we measure. Additionally, the electric and 
magnetic scattering amplitudes defined on the unit sphere with normal v 2 must also satisfy the asymptotic 
boundary conditions 


= °> where H oo = V a£ “' (13) 

We mention that this equation simply means that the electric and magnetic scattering amplitudes are 
purely tangential vector fields, possessing no normal contribution in any direction x e S 2 . 

Incoming electric and magnetic waves are ones similar to the radiating case, only they satisfy the 
conjugate of the Silver-Miiller and Sommerfeld radiation conditions. That is, rather than having constant 
phase fronts that migrate outward, in the usual sense, they travel and converge on the origin as a function 
of positive time. Additionally, by a local plane wave propagating in the direction y we mean a certain type 
of incoming wave—specifically, a Herglotz wave with a Dirac kernel supported in the direction y. 

We will now turn our attention to the case of the scattering of an incident wave pair ( E', H‘) from a 
perfect conducting scatterer. In particular, what we mean is that the tangential components of the total 
electric field must vanish throughout the surface of the scatterer, so that if v D (x) denotes this normal of the 
scatterer D at some point x, then we have 
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v D (x) a E(x) = 0, Vjc e 3D. 


(14) 


In summary, the problem is to understand how an incident pair of fields (E“, H‘) interacts with the 
conductor D and generates the scattered pair of fields (E s , H s ) as a function of various parametrizations of 
the incident wave energy, e.g., their frequency or energy and direction. We note that the problem may be 
stated and treated as a system of partial differential equations, admitting very insightful analysis and 
solution techniques via the theory of finite elements. For more on this approach see [3,4], For the purposes 
of the exposition of this report, however, we will choose to frame this problem in terms of integral 
equations for the clarity they provide in terms of the operator mappings between the scattering system 
inputs (E l , H‘) and its outputs (E s , H s ). 

In any domain Q c R within which there are no charge or current sources, each type of outgoing or 
incoming wave satisfies our frequency-dependent Maxwell system 


Va E l ’ s - ikH , s = V a H’’ ' + ikE 1 ’' = 0, k = w/c = w(|i 0 £ 0 ) *, (15) 


and hence, the total fields E and H do as well. This system of coupled partial differential equations has an 
alternative integral equation representation called the Stratton-Chu integral equation, which for a perfect 
conductor has the form 


£(•*) = E\x) + A V A j G(x-y)v(y) a H(y)ds(y), (16) 

dD 

H(x) = H‘(x) + V a |G(x-y)v(y) a H(y)ds(y). (17) 

az> 

Here, ds(y) denotes surface measure on the boundary of the conducting scatterer D, and the integrals 
exist as improper ones as x —> dD 3 y. Again, we emphasize that the pair (E , H) denotes the total 
electromagnetic and magnetic fields, and we identify the integral kernel G(x-y) appearing in both integral 
equations as the fundamental solution of the scalar Helmholtz equation. Specifically, G is the outgoing 
fundamental, or distributional, solution of 


3 

(A + k 2 )G(x-y) = 8(x-y), x,ye R 3 , A = X ( 18 ) 

j= i 
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3 

where 8 e £'(/? ) is the so-called Dirac distribution, or “delta function,” and is given by 


G(x-y) 


1 g ,fcU - yl 

4tt U~y| 


x^y. 


(19) 


The asymptotic behavior of the convolution kernel G is of extreme importance; hence, we take a few 
moments to examine it. Suppose the point y is contained within the finite ball of radius R centered at the 
origin, which we write as B R (0). Suppose further that x = x/\x\ s S 2 = {x e R 3 : x] + x 2 + x] = 1} and 
G(x) is the outgoing fundamental solution to the operator A + k 2 in R 3 . Then, 


G(x-y) 


1 e ikM 

4 n |jc| 


-ik U y) 


X —> oo . 


( 20 ) 


This claim is evident when we examine the asymptotics of the term |.t - y|. First, if y is required to 
remain in a finite region, while x may tend to infinity, then 


\*~y\ = \x\ + &£ + 0(\x\~ 


'). 


w ->■ 


( 21 ) 


Since y e B R ( 0), then 


max 

£ g gk\(x,y)\<R<oo. (22) 


Hence, 


i ik\x\ 

G(x-y)= 4^~\^ e " k(X,y \i +0(W‘)), W-»~. (23) 

The asymptotic behavior of the acoustic outgoing fundamental solution G coupled with the 
asymptotic behavior of the scattered electric and magnetic waves implies that the electric and magnetic 
scattering amplitudes are given by 


£«,(•*) = \ e ^’ ^((y) A// O0) ^xds(y) 
ao 


(24) 
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(25) 


Hjb = ~XA je lk(x,y) v((y) AE(y))ds(y). 

dD 


The formulations for the electric and magnetic scattering amplitudes now provide us an operator we 
term F, whose analytical properties will be of significant interest. It tells how and why we can image 
conducting bodies with narrowband measurements of the scattered electromagnetic waves made at a 
distance. For this reason we will examine a few more aspects of the operator F and briefly discuss how it 
indicates and describes the interaction between an incident incoming magnetic wave H‘ and a conducting 
body. More importantly, we will hit on how such a description benefits us in suggesting a reconstruction, 
or three-dimensional imaging algorithm. 

From the previous analysis we see that the electric scattering amplitude is generated by the surficial 
current j = v d a H generated on the boundary of the conductor D by the interaction of an incident pair of 
waves ( EH ‘). In summary. 


E~(x) = (Fj)(x): = J K(x,y)j(y)ds(y), (26) 

dD 

where the kernel of operator F is 


is/ * x ik —ik(x,y)* * 

K(x, y): = -—e xaxa. (27) 

471 

If we recall and write the second of the Stratton-Chu integral equations as a function of the incident 
magnetic wave 


H‘ a = J« it<y,x> a(i)ds(jc), ae7X 2 (S 2 ), (28) 

s 2 

where the subscript indicates the parametrization of H 1 —e.g., in this instance H a is a superposition of 
entire plane waves, also called Herglotz waves, and plays the role of the weighting of such square 
integrable tangential vector waves and is the Herglotz kernel—as 


H(x) = H' a (x) + V a jG(x —y)v(y) a H(y)ds(y) 

dD 


(29) 
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= : H‘(x) + (BH)(x), 


(30) 


then, formally given an known parametrization of an incident magnetic wave H ' a , the electric scattering 
amplitude is generated by 


EM) = J/fU,y)v D A((/-Br‘^)(y)^(y) (31) 

dD 


= :(F H‘ a )Cx), ieA. 


(32) 


In the next section, we present the main result of this report, which is that the coefficients in the 
generalized Fourier series representation of on A encode the boundary of the scatterer dD. Moreover, 
this result gives us a robust and efficient algorithm that can reconstruct the scatterer with limited but 
appropriately sampled values of E „ on A. The details of the properties of F and the representation of the 
signal E ^ on the array A that support these results will then be discussed in detail in Section 4. 
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3. MAIN RESULTS 


Before discussing the main results we need to address two issues. The first concerns the dependence 
of the electric scattering amplitude on the origin of the physical problem. Since we are at liberty to define 
the origin where we choose and this value changes the values of the measured quantity we have, we must 
determine this influence. Moreover, we will demonstrate that this analysis leads us to a very attractive and 
easily implementable algorithm for characterizing the scattering body. The second issue regards the vector 
representation of the tangential vector field on our measurement array. 

We begin with a quick treatment of the continuous dependence of the electric scattering amplitude 
with the origin c e R ' . In what we developed in the previous section, we assumed the origin to be the 
standard one, namely, c = 0. However, if instead we choose any arbitrary point c to be the center of our 
universe, then the electric scattering amplitude is generated by 


EM) = ( Fj(-c))Cx) 


(33) 


= J E(x, y)j{y -c)ds(y) (34) 

dD-c 

= 1 e ,k{X ’ y) xAXAj(y-c)ds(y) (35) 

dD-c 


ik ik{x, 

~4n e 


0 je ,k< ' x ’y'>xAXAj(y)ds(y) 

3D 


(36) 


= <t> c Cx)(Fj)(x), 


O c .(i) = e 


ik(x , c) 


(37) 


where dD-c is the translate of dD to the ocentered origin. Hence, given the true scattering amplitude 
produced by the physical interaction with the incident wave, we may produce its value if instead we insist 
that the point c is the origin by performing a unitary operation, i.e., a simple phase multiplication at each 
point where we measure the field. We further add that this concept of shifting and choosing various origins 
plays a very powerful and useful role in the imaging algorithm. 

We recall that the electric scattering amplitude we measure in remote sensing on a flat or curved two- 
dimensional array will consist of two orthogonal components. By this we mean that the electric scattering 
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amplitude is a vector-valued function. Also, since any flat or smoothly curved array is locally equivalent to 
the unit sphere, we begin with a very brief discussion of two-dimensional vector fields supported on the 
unit sphere S 2 . In what follows, for a scalar function w(0, <}>) defined on the unit sphere S 2 , we will use the 
so-called surficial gradient, defined by 


V 5 2 " 


1 du du \ 
sin03(|)’ dQs ’ 


(38) 


which is an element of the tangent space of the unit sphere, denoted by 7^. 

Let v be a complex-valued vector field supported on the unit sphere S 2 . Then we say that 
v e TL (S ), provided that in the representation of the vector field on the sphere S 2 


*<*> = X I 


a'y^Cx) + P'"jc A 


xeS 


n - 1 m 


(39) 


obeys the condition 


jQ<v(jc), V s2 y'"(x)> + <v(x),xa V 52 y^(x)> ^ds^x) 


< oo 


(40) 


where the terms (v, V^Y 1 ") and (v, x a V^F^ 1 ) represent the projections of the tangential field on S 2 
onto the orthonormal bases and IaV^. Here, the functions ^"(x) are the spherical 

harmonics, of order n and degree m on S 2 with x = sin0, cost)), sin0sin<|>,cos0, 0 e [0, 7 t] and 
<)) € [0, 2n ]. For our purposes we define the harmonic functions in the following way: 


4 >): = 


V 4;t (n + |m|)! 


„ 2n + 1 (n - m)\ n m, n . i 

£,„ \—\ -,-trP„(COS0)e 

m *' 4n (n + m)! n 


-l m , m > 0 
1, m < 0 


(41) 
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Above, the functions P™( cos0) are the associated Legendre functions; c.f. [5, 6] for more details on their 
various properties and governing equation. This definition of the spherical harmonics allows us to express 
their complex conjugates as simply 



Cw- 


(42) 


Also, this definition yields a complete orthonormal set on L~(S 2 ) in the sense that 


jwrfdmi) - ( 43 > 

s 2 

where 5”*’™ is the Kronecker delta. Lastly, we note that the finiteness criterion amounts to requiring that 
the generalized Fourier coefficients satisfy 


Z Z kf + |p;f <oo- (44) 

n = 1 m = -n 

We now have all the key ingredients in place and state the main result of this article as the following 
theorem. 

Theorem 1. Let B R (c) denote the sphere of finite radius R centered at the point c e R ‘ , and let E^ be 
known on some section T cz S of the unit sphere and have a series expansion in tangential spherical 
harmonics defined on T, 


E Sx) - Z Z < v /;w + P^aVs 2 );w. «r, (45) 

n = 1 m = -n 

with Fourier coefficients depending on R and c. Then the following are equivalent 

1 . B r (c) contains dD 

2. There exists a magnetic current j supported within B R (c), depending on R and c, such that 


£«(*) = J K(x-y)j(y)dy , xg T. 

Br( c ) 


(46) 
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3. The weighted energy of the Fourier coefficients, depending on R and c, is finite. 


oo n 

I I C )| 2 + |P”(R,C)| 2 . *,(R,c) = (47) 

n = 1 m = -n 

Remark 1. The theorem states several important facts all at once. In summary, each item implies one 
another, which is to say that given any item, the other two are true as well. In what follows we address 
these implications and provide some explanation. 

To begin it states that if the scatterer dD can be placed or fully contained within some sphere of radius R 
centered at the point c, then there exists a magnetic current supported on the boundary of a two- 
dimensional surface contained within this sphere, which when integrated against the kernel of the electric 
scattering amplitude generator yields the known electric scattering amplitude E x . 

Next, the result states that such a scattering amplitude must have partial Fourier coefficients that decay 
faster than (2 n/ek{R + |c|)) ”, as n tends to infinity. This is the observable quantity we exploit for various 
centers c, which then determines the value of R for each such c. Upon intersecting all such spheres for an 
appropriate collection of centers, this quantity yields an approximation of the convex hull of the scatterer, 
which if the scatterer is convex is an approximation of the scatterer itself. We note that in the limit that |c| 
tends to infinity and the approximations of the convex hull of dD become better and better, theoretically 
speaking. From a numerical point of view this approximation becomes more and more difficult, as it 
requires the computation of higher and higher Fourier modes, which is difficult due to the rapid oscillation 
of the spherical harmonics for large n. However, even numerically speaking, this approximation can still be 
done if we utilize arbitrary machine precision techniques. 

Lastly, if what we call the weighted energy of the Fourier coefficients converges for some R and center c, 
then it must be that there exists a sphere with such a radius and center that contains the totality of the 
scatterer dD, which in turn implies that the electric scattering amplitude is generated by some magnetic 
current on the surface of some two-dimensional manifold contained within such a sphere. 

In a sense, what is going on here has to do with the analytic continuation of an analytic function. The 
function is analytic in several variables, x and k, for instance. Furthermore, it is analytic in x, which 
means if we know it on some distant surface away from the scatterer dD, the function may be analytically 
continued up to the boundary of the scatterer. Where this continuation is no longer possible coincides with 
R for any given c, and means that no longer will the weighted Fourier coefficients converge, and no longer 
is there a magnetic current that can be integrated over such a smaller sphere and give rise to the known 
field 

We now offer a concise version of this result in the form of an algorithm that has been successfully 
implemented in the acoustic setting. 

Algorithm 1. (Convex 3D Imaging). 

1. Measure E ^ on some two-dimensional array A. 
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2. Determine the tangential analysis eigenfunctions spanning the array A that are used to determine 
the spectrum, i.e., generalized Fourier coefficients, of the signal E^ (functions such as V and 
x A V C in the case when A is the unit sphere 3 4 ). 

3. Optimally span a collection of points c and compute the associated generalized Fourier 
coefficients for that center. 

4. Determine the radius R for each such center by examining where the partial sums of the Fourier 
coefficients begin a rapid transition to zero. 

5. Intersect all such spheres to define the approximation of the convex hull of the scatterer. 

Remark 2. We are also compelled to make another remark concerning the uniqueness of reconstructing the 
true scatterer in this problem. Our main result states that the measured asymptotic wave we know may be 
analytically continued in a unique fashion up to some natural boundary that corresponds to the boundary of 
some scatterer. To be fair, we must mention that in general this scatterer need not coincide with the true 
scatterer of interest. Rather, what we ultimately arrive at in this process is the so-called scattering support 
of the measured field, or the scatterer. See [1,2] for more of the fundamental details concerning this object, 
as this object can be highly nontrivial and rather complex. For our purposes here, we define this object 
through the following convention. 

Definition 1. Let ch(supp f) denote the convex hull of the support of a distribution/and let s = -1/2. Then 
the convex scattering support of the far field Ef for a perfectly conducting electromagnetic scatterer is 


cS k supp E „ = p, ch(suppf), 
Ff = E_ 


(48) 


where /e 7X 2 (/?(/? ? )) 4 


What we can uniquely determine is the common element to all candidate scattering objects that 
could give rise to the measured field of interest. We need to make this interjection and state this result, 
since to date there is no rigorous proof that a perfectly conducting scatterer may be uniquely determined 


3. For more complex measurement arrays, these eigenfunctions may be analytically or numerically computed by 
treating a standard eigenvalue, eigenfunction problem whose formulation follows from the spectral theory of compact 
linear operators presented in the next section. 

4. The space of distributions ^(R 3 ) for 5 = -1/2 is the standard Sobolev space jhT ]/2 (R 3 ) of distributions having 
compact support in /?^; i.e., they are TL 2 vector fields on a two-dimensional manifold such as the surface of a sphere. 
Moreover, since we are assuming that the scatterer is perfectly conducting, we may specify s = -1/2 rather than having 
to take more general values of 5 falling within the interval (-2, 0) for which the original theory of scattering support 
was developed. 
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from one single fixed energy excitation with measurements taken on such a section of the unit sphere, or 
two-dimensional array. Nonetheless, on the basis of uniqueness results and several numerical studies of 
this problem in the acoustic setting, it is widely suspected that the convex scattering support and the 
scatterer’s true support are one and the same in the case of a perfectly conducting electromagnetic obstacle. 
A careful study and proof of this conjecture would constitute a major breakthrough in the academic realm 
of scattering theory; however, it is beyond the current scope of this report. 

Remark 3. We should also take a moment to offer a generalization of this result, which allows us to image 
more complex three-dimensional targets that need not be convex. In the case of multiply connected 
nonconvex targets, we can offer a theoretically based method that can image such systems. However, a 
practical numerical implementation of such a result has yet to be done and is beyond the work presented in 

m. 

More generally, we have the following formulation. 

Theorem 2. If Q denotes a (possibly) multiply connected finite subset of R 3 , and E m is known on some 
section T c S~ of the unit sphere and has a series expansion in tangential spherical harmonics defined on 

r, 

oo n 

S ; s r, ,49) 

n = 1 m--n 

with Fourier coefficients depending on Q then, the following are equivalent: 

1. Q. contains dD. 

2. There exists a magnetic currenty, supported within and depending on Q, such that 


£<»(*) = j’rt'C*-y)./(y)dy, xe r. (50) 

a. 

3. The weighted energy of the Fourier coefficients, now depending on Q, is finite. 


o° n 

2 >.<n> S k(Q)r + |u:(n)f. (5i) 

n - 1 m - -n 

For the purpose of brevity we shall not go into the rather complex nature of the exact asymptotics of 
these more generalized weights w n (Q), since their description is so highly connected with the geometry 
and number of components of Q. We shall have to content ourselves that this formulation extends the 
previous theorem and represents the more general methodology of imaging nonconvex scatterers, since the 
test region Q is allowed to be nonconvex, whereas the aforementioned spheres were not. 
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In the supporting theory to follow in the next section we give the appropriate indications that 
corroborate the claim made here and tie it to a result previously obtained in the scalar acoustic scattering 
case [7]. Additionally, we provide the mathematical machinery that yields a numerically feasible algorithm 
that can image such scatterers of interest. 
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4. SUPPORTING THEORY 


We begin with a basic physical treatment of the scattering problem that inspires the more rigorous 
supporting theory, which underpins the two theorems given in the previous section. 

The scattering of a fixed monochromatic electromagnetic wave of wave number k - 2n/\ from a 
perfect conductor dD situated in the ball B R (0) is equivalent to the radiation problem 


V a E s (x) - ikH\x) = V a H'\x) + ikE'\x ) = 0, x e R ? \(dD n D ), 


(52) 


where 


E *|ao = ve TL\dD) 


(53) 


is specified. The general solution of such a radiation problem in the complement of B R (0) admits the form 


E\x) = X Z °C V A ( xh { n '\kr)V:(x )) + p:v a V a (xh^\kr)Y^(x)), r = W > R, (54) 


v"'/ 


n = 1 m = -n 


where h„ is the spherical Hankel function of the first kind—c.f. [6]—and is defined by 


hn\z) = j n (z) + iy n (z) , (55) 

with j n and y n respectively denoting the standard spherical Bessel functions of the first and second kind. 

Examining the asymptotic behavior of the outgoing electric field leads to the electric scattering 
amplitude, which is known to have the form [8] 


= \ I r'-kvy^w-aysxx) 


1 - n 


(56) 
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Moreover, the tangential electric scattering amplitude is an element of TL 2 (S 2 ) only if 


X k| 2 + M 2< ~. VK>diam(dD). (57) 

n = 1 tn - -n 

The necessity of the weighted summability of the Fourier coefficients a™ and (3”' is a consequence 
of the asymptotic behavior of the spherical Hankel functions as a function of its index n. For a fixed 
positive argument z. 



n 


n —> °° . 


(58) 


Hence, has a finite L 2 norm on the surface of any sphere of radius r > R in TL 2 (B r (0 )) only if the partial 

sum of the combined Fourier coefficients decay faster than ’(Ul )| 2 as n -» °° for all r > |z|. For each 

such r, where we see this convergent behavior indicates where the radiation problem may be posed as a 
homogeneous Maxwell system, i.e., how far in toward the origin the analytic scattered field may be 
uniquely continued as a free solution to the Maxwell system. That value r* where the analytic continuation 
may no longer be performed defines the boundary of the scatterer. The analysis to follow treats the 
operators involved in the scattering problem and formalizes this basic physical principle. 



(59) 


To begin this analysis, we note that we need know only the electric scattering amplitude on some 
sector, or neighborhood, of the unit sphere. The significance of this is that there exists a collection of 
coordinate charts {<p y }, : R —» S that cover the sector of interest S 2 . What this means is that there is a 
change of coordinates so that a “piece” of the unit sphere may be made to look like a finite two- 
dimensional plane, and vice versa. Hence, in talking about the measurement array A or the unit sphere S 2 
we are talking about the same object, up to a nonlinear change of variables, so that all the norms to follow 
can be made equivalent on either the more easily treated unit sphere S or the true measurement array A. 
See chapter two of [4] for more information concerning this differential geometry. 

We proceed with a characterization of the operator that maps induced magnetic currents on the 
boundary of the scatterer to the electric scattering amplitude. 
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Proposition 1. Let dD c R 2 be smooth surface such that T (W is defined everywhere, and suppose 
there are no nontrivial solutions to the interior version of the Stratton-Chu integral equations. Then 
F : TL~(dD ) —> FL“(S“) is compact and linear and has dense range in TL 2 {S 2 ). 

Proof. Linearity is obvious as F(g, + g 2 ) = Fg x + Fg 2 for any g, 2 in TL 2 (dD). F is compact, since it has 
a so-called Hilbert-Schmidt kernel; i.e., Ke TL 2 (ds(dD) x ds(S 2 )) . Let ge TL 2 (dD ) be a magnetic 
current on dD. Let 6 and (p denote the tangential unit vectors on S 2 . Then, the previous claim follows, 

since 


11/11 


TL\ds(dD)xds(S 2 )) 




(60) 


J J 

J-2 ■'in 


167 iW'dD 


(61) 


<-^pOD)p(5 2 )<oo, 

16n 


(62) 


where p(3D) and p(S 2 ) are the finite surface measures of the scatterer and the unit sphere, respectively. 
Now, given that K has a Hilbert-Schmidt kernel, we find that given g € TL'(dD) a magnetic current on 
dD, 


2 2 2 
II ^gll 7 x 2 (5 j ) “ II ^11 TL 2 (ds(bD) x </i(S 2 ))H^ll TL 2 (dD) <0 ° ' 

2 2 

proves the image of g is included in TL (S~). 

Lastly, we note that our operator F is equivalent up to a n/2 rotation of the tangential field on S 2 

with the far-field operator presented in [8] in theorem 7.4, which is known to have dense range in 

2 2 . 7 

TL~(S ), provided that k is not a Dirichlet eigenvalue for the interior problem on dD. 

We now address the range of the operator F acting on magnetic currents on the surface of the sphere 
Br(c) in terms of the restriction of the original operator F to TL 2 current distributions supported on the 
boundary B R (c). Below, F\ dB ^ denotes this restriction. 

Proposition 2. Let be the restriction of F to B R (c). Then F^^ possesses a Hilbert 

adjoint, such that 
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= TL 2 (S 2 ) -> TL 2 (dB R (c )) 


( 64 ) 


and for g e TL 2 (S 2 ) 

( F l dB'{c)8)(y) = -( v B*(c) A v Bg(c) , yedB R (c). (65) 

Proof. Since F is compact and linear between two Hilbert spaces, then by the Riesz representation theorem 
there exists a bounded linear operator F coinciding with the action of the adjoint F* of F such that it maps 
TL (S ) to TL (dB R (c)). F may be computed by equating the pairings 

(Fji §) jl, 2 (S 2 ) ~ 8 ) TL 2 (dB K (c )) ( 66 ) 

for any j e TL 2 (dB R (c)) and g e TL~(S 2 ). Hence, we find that 

p (j>8) = ( f J’ 8\ l \ s 2 ) ( 67 ) 

= f <f e~' k{x,y) xAXAj(y)ds(y),g(x)) ids(x) (68) 

V °Br(c) 

= f <f (e~ lk{x,y) x AX Aj(y),gix))ds(x)) ids(y) (69) 

J 3B»(c) V 

= J 0B( -{ v B R (c)^\ s , e ' k(X ’ y) x AX A g(i)ds(x)J A v B r{c) ) ^ds(y) (70) 

= (j’ F 8)rL\dB R (c))- ( 71 ) 
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Hence, the claim is proved. 

Next, we address the so-called singular system, or singular value decomposition of F restricted to 
TL 2 currents on the boundary of the sphere B R (c). Since F is linear and compact then we know it admits the 
general representation 


oo n 2 

F \ 3B.(0 = Z X 0 °??«(*. c )^n!m( R > C ) ® C ) • ( 72 ) 

n = 1 ni = -n 7 = 1 

Here, the superscripts (/' = 1) and (/' = 2) refer to the tangential components in the 0 and <)> directions, 
respectively, and the term ® is the tensor or outer product of the two orthonormal bases \\i^\ n and <p 0 m for 
j = 1,2. We add that the action of F\ dB ^ act i n § on some TL 2 current j is 

oo n 2 

( F law )(i)= X X 0 c ) O', <?n! m ( R ’ • (73) 

n = 1 m - -n J ~ ^ 

We also mention that this decomposition then goes to define the Hilbert adjoint of F|*a^( C ) through 


oo n 2 

F\ *a B ,(c) = X X ® ( 74 ) 

n = 1 m = -n 7 - 1 

where the above overline denotes complex conjugation. Furthermore, we add that the countable singular 
system 


f A,<J) A>U) A,(j),°° 

Vn, m >{„ = 1, m = 




(75) 


2 2 

of F: TL~(dB R (c)) —> TL {A) any finite two-dimensional measurement array A may be constructed by 
determining those functions such that an y of the following identities are satisfied for 


j= 1,2. 


F l 


*,0> 

dB R (c)^n, m 


a 


A,(j) A,(j) 
n, m t n, m 


(76) 
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(77) 


^1 M R (c)V 


-4,0') 

n, m 


a 


A.U) A’(j) 

n, m * n, m 


171 Z7i* A, O’) / >4,0'k 2 A,(/') 

F \dB R (c) F \ *B R (c)V n X =(«„,„) V*,* 


(78) 




-4,0) 

n, m 


(a A, 0 y a,o) 

^ n, m / tr, /n • 


(79) 


Solutions of these systems may be computed numerically, with some effort, or determined through 
the computation of a collection of multidimensional integrals. These integrals stem from the coordinate 
mapping 4> from the array A and the segment of the unit sphere S 2 on which we know the scattered field 
E „; i.e., the push forward 4*: A —> S 2 is a diffeomorphism so that the pullback T -1 exists and 4> -1 : S 2 —» 
A. In this sense, the measured signal on the two-dimensional array admits the representation 


oo n 

I «:v i ,C'r' + p;v s ,v:-T-'. ( 80 ) 

n = 1 m = -n 

The next result tells us what the analysis basis for the measured electric scattering amplitude is on A; 
namely, and V also provides a more evident formulation of the eigenvalue problem, 

which must be solved to determine the singular values * and . 

Proposition 3. Let A be a smooth finite two-dimensional manifold with a boundary, such that 4' : A —» S 2 
and 4* -1 : S 2 —> A are the usual push forward and pull back mappings. Moreover, let E^ be the measured 
electric scattering amplitude on the array A and let F | 3B#(c) and F\ sb r (c) denote the restriction of F to 
dB R (c) and its Hilbert adjoint. Then 


- 4 ,( 1 ) 

VI/ 

t n, m 


= 


-1 


11/ = JC A V 

t n, m S 


y:°4>- 


(81) 


(82) 


and o 


-4.0') 

n, m 


J = 


1,2 may be determined from 
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(83) 


= (tO v s >K lo v 


-1 


^ r lae„(c)^ r l dB Ac) x 


A 


v 2 y'"°4 / = (o ' 

5 2 n ' n, 


. 4 ,( 2 ) 2 - 

_ ) X A V 
m * 




-i 


(84) 


Proof. Let E denote the electric scattering amplitude on the unit sphere S 2 and the array A, respectively, 
and let j denote the magnetic current giving rise to the two fields. Then, 


a n,m = <£t,V„ 2 C V *'"'> 


L‘(A) 


(85) 




L\A) 


( 86 ) 




(87) 


/• (1) A,(l)v 

0 > ®n, m^n, m \ 2 (R*) 


( 88 ) 


= {P x h 


o 


4,0) 'MD 

n,m t n,m 


V(/?V 


(89) 


Since Pjj is the projection of any vector field in L 2 (R\ C 3 ) 
cp^ } }, this amounts to saying 


onto the orthonormal basis 


i * A, (1) 

F \ dB R (c)Vnm 


- 4,0)„/,(l) 

®n,m tyn>m 


(90) 


holds for 


V/;°f 1 . Similarly, examining the coefficients 


in the much the same way, 


P n,,n = 


( 4 ) 


(91) 
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( 92 ) 


= (j,F\*dB R (c)i AV/T') t!(s)) (93) 


.. 4 , ( 2 ) 4 ,( 2 ). 

(j’ ^ n > m ^P n , m ) 


iV) 


/n . , 4 ,( 2 ) 4 ,( 2 ), 

<^2^ On,m Vi \ 


2 («Y 


(94) 


(95) 


implies 


F l 3B fi (c)V 


4 ,( 2 ) 
n, m 


a 


A,(2) A,(\) 
n, m t ri, m ’ 


(96) 


and, more precisely, that \|/^ 2) = jca V 5 A" ov P 1 . This completes the proof. 

The main result presented in this section owes itself in part to what is known as Picard’s Theorem. 
The theorem essentially provides a representation for a compact linear operator A between two Hilbert 
spaces H\ and H 2 in terms of the operator’s singular system, as well as a means to assess whether a given 
element of the second space H 2 is also an element of the closure of the range of A. We take a moment to 
state the theorem and refer to [8] for its proof. 

Theorem 3 (Picard). Let A : H\ —» // 2 be a compact linear operator from the Hilbert space H\ into the 
Hilbert space H 2 with the denumerable singular system {cr n , <J) n , V„}” = ,; i.e.. 


A< t>„ = o„V„ 


(97) 


and 


A A = 


(98) 


and let (.,.) denote the inner product on H 2 . Then, the equation Af = g is solvable if and only if 
* _L 

g e N(A ) and 
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oo 


< oo . 


(99) 


I 



CT 


2 

n 


Moreover, any/of the form 


f= I 

n = 1 


(8, V„) 


4> n - 


( 100 ) 


solves Af= g. 

Finally, we state the following range characterizations of F, which we will need for the proof of the 
made theorem to follow. 

Lemma 1. Let S R (c x2 ) be of two spheres of radii Rj and R 2 having centers Cj and c 2 , such that 
B R (c { ) contains B R (c 2 ) • Then 




( 101 ) 


2 

Proof. Let B R (c , 2 ) be as stated above. Let/, 2 ^ TL (dB R c^ i 2 ))- Then the trivial extension of the 
form 


j f\(x), xedB R ( Cl ) 

{ 0 , x£dB R ^c x ) 


( 102 ) 


proves the claim. 

Next, we demonstrate that in the case when the closure of B R (c,) has no intersection with the 
closure of B R (c 2 ), then the only common element in the ranges of F\^ B ) and F| ag (r) is the trivial 
electric scattering amplitude E M = 0. This amounts to saying that two disjoint conducting scatterers can 
never produce the same electric scattering amplitude on any open set away from their support. 


Proposition 4. Let B R (c, 2 ) be as before. Then 


^/?,( c i) n ^/?j( c 2) ~ ® 


) = {0} . 


(103) 
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Proof. Each tangential component of E^ on S 2 satisfies Rellich’s lemma; c.f. [8]. Hence, if we treat each 
component, this is a immediate consequence of corollary 5.1.2 in [9], 

Proof of the Main Theorem. Application of Picard’s Theorem to the identified spaces 
H\ = TL ( dB R (c )) and H 2 = TL (S ) and the range characterizations of F| 3b (c) now prove the 
theorem in the previous section. Specifically, Green’s theorem proves 1) => 2), while*the asymptotics for 
the generalized Fourier coefficients for any free solution of the Maxwell system on the complement of 
B r (c) proves 2) => 3). Lastly, assuming 3) and knowing that the range of is dense in TL 2 (S 2 ) and 

employing lemma 1 and proposition 4 we arrive at 1). (C.f. [9] for the details of these arguments.) This 
completes the proof. 

Given Picard’s Theorem, and our range characterization of the operator F| 5b (r) , we now have a 
definitive test that can determine whether a scatterer is fully within some region El of interest by means of 
testing the convergence of the sum 


ini2 


11 TL 2 (da ) 


z 


n = 1 


n 


z 


m = -n 


<£~>hO 

2 

<£~,vi 2 L> 

<C( Q > 




(104) 


In theory, if the sum does not converge, then we can conclude that the scatterer is not fully within the 
test region El, whereas if the sum does converge, we can conclude that the scattering support of the scatter 
is fully contained within this set. 


We note that this formulation states both Theorems 1 and 2 at the same time; namely, £2 may play the 
role of the sphere B R (c) or any other finite multiply connected, possibly nonconvex set. In the more general 
setting, once we have determined the appropriate singular system for the operator F restricted to the 
boundary of El —i.e., we solve the eigenvalue eigenfunction equations in proposition 3 with B R (c) simply 
replaced by the more general El —we may perform the same Picard test and conclude whether the scatterer 
is fully within El or not, as indicated by the finiteness of the series. There is good reason to believe that this 
test may be done in the electromagnetic case discussed here, again on the basis of the numerical procedure 
and results obtained in [7]. 
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5. SUMMARY AND CONCLUSIONS 


The results presented here are fundamental to many advanced radar detection and imaging efforts 
and offer us an avenue by which the three-dimensional shape of rotating bodies may be determined with 
narrowband monostatic sensors. 

We have rigorously developed and presented the underlying theory of the scattering of a single 
monochromatic, or narrowband, incident electromagnetic wave and its interaction with a perfectly 
conducting object. Moreover, we have carefully detailed how this scattering problem leads to a viable and 
robust reconstruction, or three-dimensional imaging algorithm, that is capable of estimating the size and 
shape of either a convex or multiply connected nonconvex scatterer with this extremely limited 
information. 
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